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We analyze the far field resolution of apertures which are illuminated by a point dipole located at subwavelength 
distances. It is well known that radiation emitted by a localized source can be considered a combination of 
travelling and evanescent waves, when represented by the angular spectrum method. The evanescent wave part 
of the source can be converted to propagating waves by diffraction at the aperture thereby it contributes to 
the far field detection. Therefore one can expect an increase in the resolution of objects. We present explicit 
calculations showing that the resolution at the far zone is improved by decreasing the source-aperture distance. 
We also utilize the resolution enhancement by the near field of a dipole to resolve two closely located apertures. 
The results show that without the near field (evanescent field) the apertures are not resolved whereas with the 
near field of the dipole the far zone intensity distribution shows improved resolution. This method eliminates 
the requirements of near-field techniques such as controlling and scanning closely located tip detectors. © 
2013 Optical Society of America 

OCIS codes: 050.1940 Diffraction, 330.6130 Spatial resolution, 180.4243 Near-field microscopy 



It is well known [1-3] that the use of evanescent waves 
can overcome the Rayleigh resolution limit because they 
contain high spatial frequencies and yield information 
about the smallest sub-wavelength details. However, 
evanescent waves are non-radiative and they diminish 
exponentially. They can be collected and converted to 
propagating waves by small detectors of subwavelength 
dimensions [4-6]. Such a mechanism is called near- field 
detection and requires nanomctric positioning to sam- 
ples [7] . It is analogous to tunneling microscopy by a tip 
detector. On the other hand, advances in meta-materials 
has enabled new class of lenses which compensate for the 
evanescent loss and thus restoring an image below the 
diffraction limit [8,9]. 

In the present work we propose a different approach 
to resolution enhancement. We consider an infinitesimal 
dipole located at a subwavelength distance to an aper- 
ture. Radiation from a localized source can be considered 
as a diverging spherical wave exp(ifcr)/r with suppressed 
exp(— iwt) time dependence, where k = ujjc ~ 2t:/X. 
The angular spectrum of such a wave has an infinite 
spectrum of spatial frequencies k^, ky. On propagation 
from source to the image, high-frequency components or 
evanescent part {k'^ + ^y > ^^) ^re filtered out. Only 
low frequency waves (fc^ + ky < k"^) can travel to the 
far zone. From the usual theory of diffraction, a small 
object illuminated with a propagating wave generates 
diffracted evanescent modes, and conversely, by apply- 
ing the reciprocity theorem, a small object located in an 
evanescent field converts part of this field into propagat- 
ing waves [10,11]. An aperture located in the vicinity of a 
point source diffracts the evanescent part into radiation. 
Thus, the subwavelength details of the aperture is to be 
carried to the far zone and the method do not require 



point-by-point, time-consuming scanning detection. We 
remark here that there are various physical systems to 
implement our proposal. The current technology enables 
an excited atom, molecule or a quantum dot to be used 
as a dipole [12-14]. A recent experiment [15] showed that 
a semiconductor quantum dot and a near-field coupling 
nano-antcnna could be employed to produce and detect 
single photons with almost perfect collection efficiency. 
First, we discuss the resolution of a circular aperture. 
The source is located just behind the aperture at a sub- 
wavelength distance and we consider the field at the far 
zone as depicted in Fig. 1. We assume that the aperture is 
inside a medium which is linear, isotropic, homogenous, 
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aperture (object) plane detection plane 

Fig. 1. Illustration of the coordinates for diffraction of a 
dipole field by a circular aperture. 

non-dispersive and non-magnetic. Under these assump- 
tions, all the components of the electromagnetic field are 
uncoupled and the space dependent part U{r) satisfies 
the Helmholtz wave equation at each source- free point: 



(V^ + fc2)t/(r) = 0. 



(1) 



For simplicity, we assume index of refraction equal to 1 
so that k = Lo I c. For diffraction by planar screens, equa- 
tion (1) can be solved by following Kirchhoff in applying 



Green's theorem, 



Eq. (2) becomes, 
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where the integral is taken over the plane of aperture and 
d/dn signifies a partial derivative in the outward normal 
direction at each point on S. The vectors r and R extends 
from the center of the aperture to the source and the 
screen respectively and p scans the aperture plane (see 
Fig. 1). The field amplitude U on the aperture plane 5' 
is given by the spherical wave exp(zA:|^ — p\)/\R — p\. 
If both the source and the detection plane are located 
in the far zone one can make use of the Fresnel and 
Fraunhofer approximations, \r ^ p\,\R — p\ ^ fc|p|^/2, 
and the integral given in Eq. (2) can be carried out to 
obtain, 
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where A is a constant proportional to amplitude of the 
electric field at unit distance from the source, a is the 
radius of the circular aperture, and ^ = ru/Xrz is radius 

in the spatial frequency domain. Here r\\ ~ 



the radial distance to the optical axis on the detection 
plane. The intensity distribution associated with Eq. (3) 
is referred to as the well known Airy pattern. It is rota- 
tionally symmetric around the z-axis. Figure 2(a) shows 
a cross section of the Airy pattern for an aperture size 
a = 2A. The position of the first minima, measured along 
the radial axis, is given by Ar = 0.61—^. The width Ar 
is also denoted as the Airy disk radius. The ability to re- 
solve two apertures depends on the size of the Airy disk. 
The narrower the width, Ar, is the better the resolution 
will be. For far field detection, Airy disk radius can be 
on the order of wavelength only if the size of the aper- 
ture is much larger than the wavelength itself. On the 
other hand, an aperture having a size of one wavelength 
will produce an Airy disk of radius O.Glr^, which is much 
larger than the wavelength because r^ :s> X. The effective 
size of the diffraction pattern is seen to be inversely pro- 
portional to the linear dimensions of the aperture. This is 
a direct consequence of diffraction of waves and spread- 
ing comes from the uncertainty relation Apu Afcp., > 27r. 
In far- field optics, the upper bound for Afcp , is given by 
the wavenumber k — 27r/A of the object medium because 
we discard spatial frequencies associated with evanescent 
wave components. In this case the resolution cannot be 
better than ?Tim[Ap||] = A/2 and one can reconstruct a 
low-pass filtered version of the aperture. 

Now, let us consider the setup in which the dipole is 
located very close to the aperture. In this case, Fraun- 
hofer approximations are not valid in the region between 
the dipole and the aperture. Then, the integral given in 



U{r,R) = --^ // dp,dpye'-^''-^-+''-''- 
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where $ = exp{ikrz + ik{rl+ry)/2rz} is a phase factor. 
Our aim is to compare this result with Eq. (3). To ana- 
lyze the angular spectrum of the field one can expand the 
spherical wave into an angular spectrum of plane waves 
as given by Weyl's [16] integral. 
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fc2 -fc2 for k?,>k^. 



(5) 
(6) 



Rz = — e is the source-aperture distance on the z-axis 
and /cm = k^ + ky. By substitution from Eq. (6) and 
Eq. (5) into Eq. (4) it can be verified that the angular 
spectrum of U{r, R) has the following transverse wave 
vector component 



= fc- 



/C||. 



(7) 



This contains the source wave vector multiplied by the 
numerical aperture (r^\/rz) and the lateral wave vector, 
fc||, of the dipole. Therefore the spatial frequency exceeds 
that of Eq. (3) by /c] | . In principle, there is no upper limit 
for ^11 because a point source in space is characterized 
by an infinite spectrum of spatial frequencies kx, ky. A 
spherical wave is well localized in position but its mo- 
mentum has a direction that is totally uncertain. For 
large e, high frequency components are filtered out and 
we get, K|| < fc, i.e. the usual diffraction limit. On the 
other hand, for e < A, the contribution coming from the 
higher spatial frequency modes, {k?, > k^) can not be ig- 
nored. It is possible to collect higher spatial frequencies 
partially by locating the dipole closer to the aperture. 
Note that even for k?. = k^ + ky > fc^, the wave in Eq. 
(4) is propagating because of the phase factor $. 

To see the effect of the near field of the dipole we 
numerically calculate the intensity associated with the 
integral given in Eq. (4). Figures 2(a)-2(d) represent the 
intensity of diffracted field for a given radius of the aper- 
ture and for different values of dipole-apcrture distance, 
e. The bandwidth of spatial frequencies increases as e de- 
creases, because the contribution coming from larger fcu 
values increases as e gets smaller in the integral given by 
Eq. (5). Comparing Fig. 2(a) to Figs. 2(5) — 2(d) we can 
understand how the near field of the dipole increases fc| 
and changes the resolution of the object. Higher spatial 
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Fig. 2. (Color online) Far field image of a circular aper- 
ture of radius a ~ 2A. In (a) the dipole is located in the 
far zone and the intensity pattern is an Airy disk (solid 
line) . To check the validity of Fraunhofer approximations 
we also show the result (dashed line) of Kirchoff integral 
given in Eq. (2). Figures (6) to (d) show intensities with 
decreasing dipole-aperture distance, e. rz is lOOA in all 
figures. 

frequencies, fcn > fc^, which are totally lost in classical 
imaging (see Fig. 2(a)) are partially converted into prop- 
agating waves. What we are showing is in a sense inverse 
of what one does in Fresnel diffraction as we keep the 
source distance much smaller than a wavelength. How- 
ever the observation plane is in far field. The precise 
mathematical connection between our results and fres- 
nel diffraction needs to be examined further using the 
rigorous wave equations. 

The far field resolution enhancement provided by the 
evanescent component of the dipole can also be consid- 
ered for the diffraction by two closely located apertures. 
We calculate the far field using the same scalar diffrac- 
tion integral in Eq. (4) . This time the integration surface 
is chosen to be two circular holes as illustrated in Fig. 3. 
We choose the radii of the apertures to be A and A/2. We 
consider two different separation distances, d = 2A and 
d = A/4, measured from the rim. First we calculate the 
patterns for the source located in the far zone for com- 
parison (see Figs. 3 (a) and 3 (d)). The images look like 
an interference pattern from two slits and the objects are 
not resolved. As the source gets closer in the subwave- 
length range, the objects appear as two separated spots 
and they are resolved (see Figs. 3 (5), (c) and 3 (e), (/)). 

To summarize, we have described how the high spatial 
frequencies of a localized scalar field can be utilized to 
enhance the resolution of objects in the far zone. We use 
the scalar diffraction theory and the method of angu- 
lar spectrum representation to compute the images nu- 
merically. The principle of enhanced image resolution is 
explained by analyzing the transverse wave vector com- 
ponent of the the diffracted field in the object plane. 
We also demonstrate how partially collected evanescent 
waves by the scalar dipole can resolve two closely lo- 
cated objects. The principle of resolution enhancement 
by using evanescent waves in this paper do not require 
a scanning probe or near-field detection schemes. There- 
fore, the far-field microscopy techniques can be employed 



for the image reconstruction. 
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Fig. 3. Computed far field images of two separated cir- 
cular apertures, r^ is 500A in all figures. The insets in 
(d), (e) and (/) show the intensities in two dimensions. 
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